HD-THEP-06-07 



Stability and Symmetry Breaking in the General 
Two-Higgs-Doublet Model 
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A method is presented for the analysis of the scalar potential in the general Two-Higgs-Doublet 
Model. This allows us to give the conditions for the stability of the potential and for electroweak 
symmetry breaking in this model in a very concise way. These results are then applied to two 
different Higgs potentials in the literature, namely the MSSM and the Two-Higgs-Doublet potential 
proposed by Gunion et al. The known results for these models follow easily as special cases from 
the general results. In particular, for the potential of Gunion et al. we can clarify the stability and 
symmetry breaking properties of the model with our method. 
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The Standard Model (SM) of particle physics is the- 
oretically consistent and experimentally successful to 
date The recently observed neutrino masses are 

very small and can be neglected in most high-energy ex- 
periments. Only one ingredient of the SM, the Higgs 
boson, has yet to be discovered. From direct searches 
at LEP a lower bound on the Higgs-boson mass of 
114.4 GeV at 95% C.L. is obtained when the data from 
the four LEP collaborations are combined Fur- 
thermore, from measurements of electroweak precision 
observables at LEP, SLC and NuTcV that depend on 
the Higgs-boson mass through radiative corrections and 
from other H^-boson measurements, one obtains (see 
Tab. 10.2 in l\) the prediction m H = Qlttl GeV - Since 
the one-loop corrections depend on the Higgs mass only 
as log(mjy /mw) the errors are rather large and the up- 
per error is larger than the lower one. Such an indirect 
determination of a particle mass is very successful in case 
of the top-quark, see also Tab. 10.2 of [lj. However, there 
the observables have a quadratic dependence on the mass 
and are therefore much more predictive. A direct discov- 
ery of the Higgs boson at the LHC is presumably possible 
up to a mass of about 1 TeV, see e.g. One or several 
Higgs bosons may be found opening up the direct study 
of the scalar sector of particle physics. 

Despite its experimental success, the SM is not sat- 
isfactory as a fundamental theory not only because it 
contains a large number of parameters, that is coupling 
parameters and particle masses: The squared physical, 
renormalised Higgs-boson mass m 2 H , which we expect 
to be of the order of the squared vacuum expectation 
value v 2 (250 GeV) 2 of the Higgs field, receives large 
quantum corrections. These corrections depend quadrat- 
ically on the particle masses that the Higgs boson cou- 
ples to. This means that the Higgs-boson mass is sensi- 
tive to the heaviest particles of the theory, for instance 
from physics which alter the high-energy behaviour at 
the GUT or the Planck scale. In principle, the corrections 
could be much larger than v 2 but cancel with the squared 
bare mass so that the difference gives m 2 H ~ v 2 . However, 
such a fine-tuning is usually considered unnatural. For 
this so-called naturalness problem see e.g. 5]. A system- 
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atic cancellation of quantum corrections to the squared 
Higgs-boson mass is provided by supersymmetry |(j . The 
simplest supersymmetric extension of the SM is the Mini- 
mal Supersymmetric Standard Model (MSSM) 0, which 
has been studied extensively in the literature. In the 
MSSM one has two Higgs doublets. In further exten- 
sions like the Next-to-Minimal Supersymmetric Standard 
Model (NMSSM) a further Higgs singlet is added. 

In this paper we study a class of general models hav- 
ing a scalar sector with two Higgs doublets. We sup- 
pose that the SU(3)c x SU(2) L x U(l)y gauge symme- 
try holds. In its simplest version the fermion content of 
such a model is assumed to be the same as that of the SM. 
The same is assumed for the gauge bosons, thereby avoid- 
ing to introduce new fundamental interactions. In prin- 
ciple, electroweak symmetry breaking (EWSB) works in 
these models in a similar way as in the SM. The La- 
grangian contains terms that consist only of scalar fields 
without derivatives. These terms form the scalar po- 
tential at tree-level and are responsible for the stability 
and symmetry-breaking pattern of the model. Further, 
through their covariant derivatives the scalars have cou- 
plings to the gauge bosons, and through Yukawa interac- 
tions they couple to fermions. After EWSB these terms 
are responsible for the generation of the gauge-boson and 
fermion masses, respectively. With increasing number of 
scalar fields the number of parameters in the potential 
becomes soon very large. For instance, as we shall see 
in Sect. |3| below, there are 14 parameters to describe the 
most general potential with two Higgs doublets in con- 
trast to only two parameters for one doublet. Therefore 
the characterisation of the symmetry breaking for dif- 
ferent regions in parameter space becomes increasingly 
complicated. We present a formalism for the analysis of 
stability and of spontaneous symmetry breaking in mod- 
els with two Higgs doublets. 

There exists a vast amount of literature on the Two- 
Higgs-Doublet Model (THDM) , where typically the num- 
ber of parameters of the potential is restricted by contin- 
uous or discrete symmetries. For instance in a detailed 
discussion of the symmetry-breaking pattern for different 
regions in parameter space is given for the THDM where 
a symmetry is imposed on the Higgs potential and 
in the stability of the CP conserving THDM is in- 
vestigated. We also want to mention an approach 11] to 
deduce the parameter constraints from stability and sym- 
metry breaking in one specific model, the THDM intro- 
duced by Gunion et al. ( 1131 . see chapter 4). Further, in 
two complementary works |l3l [l4| , the hierarchy between 
the charge breaking and the charge conserving minima is 
investigated for the general THDM. Basis independent 
techniques for t he g eneral THDM are used in different 
recent works [ll^ll6tlr^ll8| to analyse various aspects of 
the vacuum. 

Here we deduce the parameter constraints from the sta- 
bility and from the electroweak symmetry breaking con- 
ditions in the general THDM. The global minimum of the 
potential is found by the determination of all stationary 



points. Our results agree with those of 9] if we impose 
the conditions on our parameters such that the potential 
is invariant under that discrete symmetry. Moreover, our 
formulation of the criteria for stability and EWSB of the 
potential is very concise and should, therefore, be inter- 
esting for its method. This general method, where the 
potential is expressed in terms of gauge invariant func- 
tions, was proposed already in a previous work [19|. 

We also remark that the scope of the present analysis 
is the classical level. In a more detailed study quantum 
corrections should be taken into account. Some aspects 
of radiative corrections for the Higgs potential in con- 
strained n-Higgs-Doublet Models are discussed in [20| . 
The question if stability at the classical level is really 
necessary for a consistent quantum theory was put for- 
ward a long time ago by Symanzik |21| . The answer given 
in |22| was that, indeed, it is necessary. Thus, the results 
obtained at the classical level are important for the full 
theory. 

This work is organised as follows: In Sect. [21 we give 
general motivations for an extended scalar sector and re- 
view theoretical and experimental constraints on Higgs- 
boson masses in models with two Higgs doublets. In 
Sect. |3] we present the Lagrangian for the THDM. We 
introduce our notation for the Higgs potential, which is 
expressed in terms of gauge invariant functions of the 
fields. In Sect. 01 we analyse the conditions for the stabil- 
ity of the potential. In Sect. we derive expressions for 
the location of the stationary points of the potential. The 
conditions derived from spontaneous symmetry breaking 
of the electroweak gauge group SU(2)l x U{l)y down 
to the electromagnetic gauge group U(l) em are given in 
Sect. In Sect. we specify the potential after EWSB 
in our notation. Eventually, in Sect.|H|the results are ap- 
plied to two specific models with two Higgs doublets, the 
MSSM and the model of Gunion et al [lj- We present 
our conclusions in Sect. [5J In Appendices 1X1 and FBI we 
discuss the structure of the space of gauge orbits for the 
THDM and for the general model with arbitrary number 
of Higgs doublets, respectively. 



2. MOTIVATIONS FOR AN EXTENDED 
HIGGS SECTOR 

Given the fact that theoretically the mechanism 
of EWSB in the SM with one Higgs doublet is working 
well and that experimentally not even one fundamental 
scalar particle is discovered yet, what are the motivations 
to consider an extended Higgs sector? Some reasons are 
as follows. 

A promising candidate for a theory that solves the 
naturalness problem and has a higher symm etry than 
the SM is the MSSM 0, for reviews see e.g. [HHl. Also 
the MSSM contains fundamental Higgs fields that are re- 
sponsible for the generation of masses. As the minimal 
supersymmetric extension of the SM, the MSSM has two 
scalar Higgs doublets, being the minimum for an analytic 
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superpotential and the absence of triangle anomalies. An 
extended Higgs sector can improve gauge-coupling unifi- 
cation at high scales (2]|- In particular, supersymmetric 
models allow the unification to occur at a sufficiently high 
scale consistent with the non-observation of proton de- 
cays [2(t| . We remark that supersymmetry imposes many 
relations between the parameters of the potential of the 
most general model with two doublets. Cosmology pro- 
vides an additional reason for a non-minimal Higgs sec- 
tor. The experimental lower bound on the Higgs-boson 
mass in the SM, m# > 114.4 GeV, is too high for the 
electroweak phase transition in the early universe to pro- 
vide the thermal instability that is necessary for baryo- 
genesis (21), for a review see [28|. In this respect models 
with additional scalar particles are more promising than 
the SM, see |2jj. Last but not least, given the large spec- 
trum of fermion masses and the fact that the fermion- 
scalar interactions are responsible for their generation, 
the idea does not seem too abstruse that several scalar 
particles are involved in this mechanism. There are three 
known generations of fermions so why should there exist 
only one Higgs boson? 

A study of the general Higgs sector of a theory pos- 
sessing the gauge group SU(3)c x SU(2)l x U(1)y was 
presented in |29j . In the following the strong interaction 
gauge group SU(3)c will play no role, so we shall not 
mention it further. In [2^| the scalar fields, collectively 
denoted by Xj are supposed to transform under a gen- 
eral representation of the gauge group SU(2)l x U(1)y ■ 
Such a representation may be reducible and consists of 
complex unitary and real orthogonal parts. However one 
can show that without loss of generality it can be as- 
sumed that x carries a real orthogonal representation of 
SU(2) L x U{l) Y 0. For the THDM this correspon- 
dence is demonstrated in Appendix B of [3(J • The scalar 
potential is then assumed in [2|j to have a non-zero vac- 
uum expectation value 



^ = (0|x|0> ^ 



(2.1) 



and to leave the electromagnetic subgroup t/(l) em un- 
broken as usual. We use the boldface letter here in order 
to signify that v is, like in general a multi-component 
vector. One can then compute particle masses and cou- 
plings for arbitrary representations of scalars. However, 
only some representations are allowed in order to be in 
agreement with experimental data. 

One main restriction originates from the observed high 
suppression of flavour-changing neutral currents. A way 
to ensure this in the theory is to require that all quarks 
of a given charge receive their masses from the vacuum 
expectation value of the same Higgs boson |3l| ■ Since we 
analyse only the scalar potential in this work and do not 
specify the Yukawa interactions we shall not discuss this 
condition further here. 

Very relevant for us here are the consequences for 
the Higgs sector obtained from the accurately measured 
p-parameter, which relates the masses of the W and 
Z bosons, mw and mz, to the weak mixing angle 6 W : 



mw 



cos 6 W mz 



(2.2) 



Experimentally the p-parameter is very close to one |32j | 
and this suggests to require theoretically p = 1 at tree 
level. This is indeed the case for the SM. For the most 
general Higgs model as studied in |2^| one finds the fol- 
lowing. 

It is convenient to extend the real representation car- 
ried by the general Higgs field ^ to a unitary representa- 
tion of the same (complex) dimension and to decompose 
it into representations with definite values (t,y), where t 
and y are the weak-isospin and weak-hypercharge quan- 
tum numbers, respectively. We have 



t 



(2.3) 



and, for reasons discussed in [23, suppose the hyper- 
charge quantum numbers y to be rational numbers. The 
normalisation is such that the charge, hypercharge and 
third component of weak isospin matrices are related by 



Q = T 3 + Y. 



(2.4) 



Then the squared gauge-boson masses, (see (2.43) 
in [2jj ), are given by 



' w 



sm( 



sin #w cos 6 



Y, [^ + i)-y 2 ] v T nt, y)v, 

t,y 

(2.5) 

2 

Y,V 2 v T P{t,y)v, (2.6) 

t,y 



where P(i, y) is the projector on the subspace with rep- 
resentation (t,y). Here, the positron charge e, and the 
sine and cosine of the weak mixing angle are defined in 
terms of the gauge couplings g and g' as in the SM, see 
for instance 1331: 



sin ( 



cosf 



9 



vV 



9' 



fj 



Vg^+g 7 ^ 

gs'm8 W7 



where we use the same notation as in [2£ 
in HI that 



v T I>(t,y)v ^ 



is only possible if 



ye {-*,-*+ !,...,*}. 



(2.7) 
(2.8) 

(2.9) 
It is shown 

(2.10) 
(2.11) 
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Inserting the expressions for mw and mz 
tion 12. 2[) one obtains [2!j 



P = 



Et, y [t(t + 1) - V 2 } v T TP(t,y)v 



Et, v 2 y 2 v T V(t,y)v 



in the defini- 



(2.12) 



To obtain p = 1 one can either fine-tune the parameters 
of the potential in order to get the right vacuum expecta- 
tion values, which seems rather unnatural and is therefore 
not considered here. Or one can only allow those repre- 
sentations in (|2.5|l and l|2.6[) that separately lead to p = 1. 
There are infinitely many such representations 34] , start- 
ing with the doublet with t = 1/2 and y = ±1/2, and the 
septuplet with t = 3 and y — ±2. From each of these rep- 
resentations one or more copies are allowed and one still 
gets p = 1. Furthermore, the singlet with y — and all 
representations with 



yi {-t,-t+ !,...,*} 



(2.13) 



can occur because they do not contribute to the sums 
in 

The simplest possibility to extend the Higgs sector of 
the SM keeping p = 1 at tree level is, therefore, to allow 
for more than one Higgs doublet. In these models the 
shape of the scalar potential depends on many parame- 
ters and can be quite complicated. As mentioned above 
it is the potential that is responsible for the scalar self- 
interactions and — together with the interaction terms of 
the scalars with the respective particle — for the gener- 
ation of the masses. Therefore one is interested in the 
conditions that one has to impose on these parameters 
in order to render the potential stable and to guarantee 
spontaneous symmetry breaking from SU(2) L x U{l)y 
to £7(1 ) cm . In this paper we consider two Higgs dou- 
blets, that is the THDM. We remark that after EWSB 
three degrees of freedom in the scalar sector reappear 
as longitudinal modes of the massive gauge-bosons. All 
other degrees of freedom of the scalar sector correspond 
to physical Higgs bosons, that is with each additional 
doublet four (real) physical scalar degrees of freedom are 
added to the model. In the THDM, there are altogether 
five physical Higgs particles: Three neutral Higgs bosons 
h°, H° (where conventionally m^o < m^a) and A , as 
well as two charged Higgs bosons . If the Higgs po- 
tential is CP conserving the neutral mass eigenstates are 
also CP eigenstates, where h° and H° are scalar bosons 
and A is a pseudoscalar. There exist various studies of 
the phenomenology of the THDM in the literature, for 
an overview and further references see for instance |35j . 

For the MSSM a large number of Feynman rules involv- 
ing Higgs bosons is derived in j3||. The phenomenology 
of the Higgs bosons in the MSSM is further developed 
in [33|. We remark that in models that possess an ex- 
tended Higgs sector (and may also contain further non- 
SM particles) for certain regions of the parameter space 
there often exists one neutral Higgs boson that behaves 
similarly to the SM Higgs boson. For instance, the MSSM 
Higgs sector is described by two parameters, which can 



be chosen as the mass of the pseudoscalar boson m A o and 
the ratio tan /3 of the vacuum expectation values of the 
two Higgs doublets. In the decoupling limit m A o 3> mz, 
where practically m A o > 200 GeV is sufficient [3^, one 
neutral Higgs boson h° is light and has the same cou- 
plings as the SM Higgs boson whereas the other Higgs 
bosons H°, A and i? ± are heavy and decouple. If there 
exist light supersymmetric particles that couple to h it 
may be comparatively easy to distinguish h° from the 
SM Higgs boson even if it has SM like couplings; this 
is because h° can decay into the light supersymmetric 
particles if kinematically allowed. Further, if the light 
supersymmetric particles couple to photons (gluons) the 
one-loop 77/1 (ggh°) coupling is modified by their con- 
tribution to the loop, thus the branching ratios of h° 
differ from those of the SM Higgs boson. If all heavy 
Higgses are beyond kinematical reach in the decoupling 
limit, such precision measurements are the only way to 
distinguish h Q from the SM Higgs boson. Notice that at 
an e + e~ collider like the ILC [33 the heavy Higgs states 
can only be produced pairwise in the decoupling limit so 
that the kinematical limit may be very crucial. However, 
at a 77 collider s-channel resonant H and ^4° produc- 
tion |4(| is possible so that only the available cm. energy 
of the 77 system limits the masses which can be explored. 
In the 77 option of an ILC one expects that the maximal 
useful 77 cm. energy will be 80% of the cm. energy in 
the e + e~ mode p4| . 

Present experiments give the following exclusion re- 
gions for various versions of the THDM. The OPAL col- 
laboration has performed a parameter scan for the CP 
conserving THDM and excluded at 95% C.L. large 
parts of the region where 



1 GeV < m h0 < 130 GeV, 
3 GeV < m A o < 2 TeV, 
0.4 < tan/3 < 40, 

7T 7T 7T 7T 
a — , , 0, — , — . 

2' 4' ' 4' 2 



(2.14) 



Here a is a mixing angle for the two states h° and H°. 
Further, the approximate region where 



1 GeV < m h o < 55 GeV, 
3 GeV < m A o < 63 GeV 



(2.15) 



is excluded for all tan/3 values for negative a. In a com- 
bined analysis 0] of the four LEP collaborations a lower 
bound on the mass of the charged Higgs in models with 
two Higgs doublets like the THDM or the MSSM, ap- 
proximately 



m H ± > 78.6 GeV 



(2.16) 



is determined. In another analysis |4J| of the four LEP 
collaborations signals for neutral Higgs bosons at differ- 
ent benchmark points of the MSSM were searched for. 
Here the limits 



m h o > 91.0 GeV, 
m A o > 91.9 GeV 



(2.17) 
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at 95% C.L. are obtained. Under the assumption that 
the "left-right" -stop mixing is maximal and with con- 
servative choices for other MSSM parameters the region 
0.5 < tan/3 < 2.4 is excluded at 95% C.L. 



3. THE GENERAL TWO-HIGGS-DOUBLET 
MODEL 

We denote the two complex Higgs-doublet fields by 

(3-D 

with i — 1,2. Hence we have eight real scalar degrees 
of freedom. The most general SU(2)l x U(1)y invariant 
Lagrangian for the THDM can be written as 

i^THDM = + -^Yuk + -Sf', (3.2) 

where the pure Higgs-boson Lagrangian is given by 

^ V =J2 V>»Vi)Hv*<Pi)-V{<px,<p2). (3.3) 

i=l,2 

This term replaces the kinetic terms of the Higgs boson 
and the Higgs potential in the SM Lagrangian. The co- 
variant derivative is 

XV = M + wW^T a + ig'B^Y, (3.4) 

where T a and Y are the generating operators of weak- 
isospin and weak-hypercharge transformations. For the 
Higgs doublets we have T a = r a /2, where r a (a — 1,2,3) 
are the Pauli matrices. We assume both doublets to 
have weak hypercharge y = 1/2. Further, ~S?Yuk are 
the Yukawa-interaction terms of the Higgs fields with 
fermions. Finally, _Sf' contains the terms of the La- 
grangian without Higgs fields. We do not specify -Sfyuk 
and Jzf' here since they are not relevant for our analysis. 
The Higgs potential V in the THDM will be specified 
below and discussed extensively. 

We remark that in the MSSM the two Higgs doublets 
Hi and H 2 carry hypercharges y = —1/2 and y = +1/2, 
respectively, whereas here we use the conventional defini- 
tion of the THDM with both doublets carrying y = +1/2. 
However, our analysis can be translated to the other case, 
see for example (3.1) in [3(|, by setting 



where e is given by 

(3. 6 ) 

The most general gauge invariant and renormalisable 
potential VQpi, (p 2 ) for the two Higgs doublets (fx and (p 2 



is a hermitian linear combination of the following terms: 

ViVj* {vW^iVkVl), (3-7) 

where i, j, k, I € {1, 2}. It is convenient to discuss the 
properties of the potential such as its stability and its 
spontaneous symmetry breaking in terms of gauge in- 
variant expressions. For this purpose we arrange the 
SU(2)l x U(l)y invariant scalar products into the her- 
mitian 2x2 matrix 

K:= (^f 1 ^fi) (3 . 8) 

and consider its decomposition 

£y = i (KoSij+Kaafj), (3.9) 

using the completeness of the Pauli matrices a a (a = 
1,2,3), together with the unit matrix. The four real co- 
efficients defined by the decomposition (|3.9|) are given by 

K = l pli Pi , tf„ = fataK-, (a = 1,2,3). (3.10) 

Here and in the following summation over repeated in- 
dices is understood. Using the inversion of H3.10JI . 

<p\<px = (K + K 3 )/2, tp\<p 2 = (K t + iK 2 )/2, 

<pl<pi = {K - K 3 )/2, y\<px = [K x - iK 2 )/2, 

(3-H) 

the most general potential can be written as follows: 

V(<p 1 ,<p 2 ) = V 2 +V i , (3.12a) 
V 2 =toK + taK ai (3.12b) 
Vi = VooK% + 2K oVa K a + K aVab K b , (3.12c) 

where the 14 independent parameters £oi £a, ?7oo> Va and 
Vab = Vba are real. We subsequently write K := (K a ), 
€ ■= (60. V ■= (Va) and E := (77^). 

Now we consider a change of basis of the Higgs fields, 
<pi — > tp^, where 

($-(££:)(£)■ <»■»> 

Here 

is a 2 x 2 unitary transformation. With (|3.13() the gauge 
invariant functions (|3.10|l transform as 

K' = K , K' a = R ab (U)K b , (3.15) 

where R a b(U) is defined by 

U ] a a U = R ab {U) a b . (3.16) 
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The matrix R(U) has the properties 

R*(U)=R{U), R T {U)R(U) = 1, detR(U) = l, 

(3.17) 

where 1 denotes the 3x3 unit matrix. That is, R(U) £ 
SO (3). The form of the Higgs potential (|3.12|) remains 
unchanged under the replacement 13.15fl if we perform 
an appropriate transformation of the parameters 

rim = mo, v' = R(U)v, (3.18) 

E' = R(U)ER T {U). 

Moreover, for every matrix R with the properties l|3.17|l , 
there is a unitary transformation l|3.13|l . We can there- 
fore diagonalise E, thereby reducing the number of pa- 
rameters of V by three. The Higgs potential is then de- 
termined by only 11 real parameters. 

The matrix K_ is positive semi-definite, which follows 
immediately from its definition Ij3.8|l . With Kq — tr K 
and Kq — K 2 — 4 det K_ this implies 

K > 0, K 2 - K 2 > 0. (3.19) 

On the other hand, for any given Kq, K fulfilling l|3.19|l . 
it is possible to find fields (fi obeying (|3.1U|) . Further- 
more, all fields obeying (|3.10|l for a given Kq,K form 
one gauge orbit. This is shown explicitly in Appendix El 

Thus, the functions Ko,K a parametrise the gauge or- 
bits and not a unique Higgs-field configuration. Specify- 
ing the domain of the functions Kq, K a corresponding to 
the gauge orbits allows to discuss the potential directly 
in the form (|3.12[1 with all gauge degrees of freedom elim- 
inated. It is curious to note that the gauge orbits of the 
Higgs fields of the THDM are parametrised by Minkowski 
type four vectors (Kq, K) which have to lie on or inside 
the forward light cone. 

In the following sections we derive bounds on the pa- 
rameters of the potential that result from the conditions 
that 

• the potential V is stable, 

• we have spontaneous symmetry breaking of 
SU(2) L x U(1) Y down to U(l) em . 



4. STABILITY 

According to Sect. we can analyse the properties of 
the potential (|3.12|) as a function of Kq and K on the do- 
main determined by K > and Kq > K 2 . For Kq > 
we define 



k := K/Kq 



(4.1) 



In fact, we have Kq = only for tpi = ip 2 = 0, and the 
potential is V = in this case. From (|3.12|) and l|4.1|) we 
obtain for Kq > 



V 2 = K J 2 (k), J 2 {k) 
Vi = K 2 Ji{k), J 4 (k) 



£o+£ T fc, (4.2) 
r/oo + 2rj T k + k T Ek, (4.3) 



where we introduce the functions J 2 (k) and Ji(k) on the 
domain \k\ < 1. 

For the potential to be stable, it must be bounded from 
below. The stability is determined by the behaviour of V 
in the limit Kq — > oo, hence by the signs of J^k) and 
J 2 (k) in l|4.3(l and (|4.2jl . For the theory to be at least 
marginally stable 



J 4 (fc) > or 

J 4 (k) = and J 2 {k) > 



for all \k\ < 1 (4.4) 



is necessary and sufficient, since this condition is equiv- 
alent to V > for Kq — > oo in all possible directions k. 
The more robust stability property V — > oo for Kq — > oo 
and any k can either be guaranteed by 



J 4 (fc) > or 



J 4 (fe) = and J 2 (k) > 
in a weak sense, or by 

J 4 (k) > for all |fe| < 1 



for all |fc| < 1 (4.5) 



(4.6) 



in a strong sense, that is by the quartic terms of V solely. 

To assure J 4 (fe) is positive (semi-)definite, it is suf- 
ficient to consider its value for all stationary points of 
J 4 (fc) on the domain |fe| < 1, and for all stationary points 
on the boundary |fe| = 1. This holds, because the global 
minimum of the continuous function J 4 (fc) is reached on 
the compact domain |fc| < 1, and it is among those sta- 
tionary points. This leads to bounds on 7700, r\ a and rjaj,, 
which parametrise the quartic term V 4 of the potential. 
For |fc| < 1 the stationary points — if there are any — must 
fulfil 

Ek = -T] with |fe| < 1. (4.7) 
If det E 7^ we explicitly obtain 

Mk)\ stat =Voo-V T E- 1 V if 1-V T E- 2 V>0, (4.8) 

where the inequality follows from the condition |fc| < 1. 
If det E = there can exist one or more "exceptional" 
solutions k of (|4.7|l . They, again, have to obey |fe| < 1. 
For |fe| = 1 we must find the stationary points of the 
function 

F 4 {k,u) := J 4 (fc)+7i(l-fc 2 ), (4.9) 

where u is a Lagrange multiplier. Those are given by 

{E-u)k=-ri with |fc| = 1. (4.10) 

For regular values of u such that det(E — u) ^ the sta- 
tionary points are given by 



(4.11) 



and the Lagrange multiplier is determined from the con- 
dition k T k = 1 after inserting (|4. 1 1|1 : 



1 - T] T (E- U y 2 ri = 0. 



(4.12) 
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We thus obtain the solution 



M k )\ st at =U + T]Q -T] T (E-U) 1 T], 



(4.13) 



where u is a solution of l|4.12Jl . Also for |fc| = 1, de- 
pending on the parameters r\ a and r\ a bi there can be ex- 
ceptional solutions (fc, u) of l|4.10|l where det(E — u) = 0, 
i.e. where u is an eigenvalue of E. 

The regular solutions for the two cases |fc| < 1 and 
|fc| = 1 can be described using one function only. Con- 
sidering (|4.9(l and l|4.11|l . we define 



f(u) := F 4 (k(u),u), 

with k(u) as in (US) . This leads to 

/(it) = u + rjoo - V T ( E - u)~ 
f(u) = l- V T (E-u)- 2 r 1 , 



(4.14) 



V, 



(4.15) 
(4.16) 



so that for all "regular" stationary points k of Ji(k) 



/(«) 
/'(«) 



Mk)\ t 
i-k 2 



(4.17) 
(4.18) 



holds, where we set it = for the solution with \k\ < 1. 
There are stationary points of J^(k) with |fc| < 1 and 
|fc| = 1 exactly if /'(0) > and f'(u) = 0, respectively, 
and the value of J4(fe) is then given by f{u). 
In a basis where E = diag(/Ui, /12, ^3) we obtain: 



f(u)=u + r, 0Q -J2- 



a— 1 



/'(«) = i-E 



^ (Ma - W) 2 



(4.19) 
(4.20) 



The derivative f'(u) has at most 6 zeros. The shape 
of f(u) and f'(u) for a (purely didactical) set of parame- 
ters where f'(u) has 6 zeros can be seen in Fig.^ Notice 
that there are no exceptional solutions if in this basis all 
three components of rj are different from zero. 

The function f(u) given by l|4.15J) allows us to discuss 
also the exceptional solutions of (|4.7fl and (|4.1(J|) . Con- 
sider first |fe| < 1 and suppose that detE — 0. In the 
basis where E is diagonal we have 



and l|4.7(l reads 



det E = fii [12 M3 = 



j«i fci = -rji, 

M2 k 2 = -T)2, 
M3 h = -»73- 



(4.21) 



(4.22) 



Clearly, a solution of l|4.22|l is only possible if with \i a = 
also r\ a — (a = 1,2,3). Therefore, we see from (I4.19|) 
that exceptional solutions with |fe| < 1 are only possible if 
f(u) stays finite at u = 0. That is, the pole which would 
correspond to /i a = must have residue zero. Suppose 



0.5 - 




-0.5 - 



0.01 0.02 



0.03 0.04 




0.01 0.02 0.03 0.04 



FIG. 1: The stability determining functions 

f'(u) and f(u) as given by 14.201 and (14.191 
with 7/oo=0.05, (m, fi 2 , (is) = (0.01,0.02,0.03) and 
772,773) = (0.002,0.002,0.002). 



now that indeed rj a — for all a where fi a = 0. Take as 
an example = [1% = and 771 = 772 = 0, but 113 7^ 0. 
Then we get the general solution of l|4.22|l as 

, V3 
K3 = , 

M3 

with k\, k 2 arbitrary but satisfying 

2 



(4.23) 



k 2 



We can write this as 



< 1. 



fe = fen + fc_|_, 



where 



k \\ = — V, 

M3 



Ek ± = 0, 

2 



fci < 1 - fcf = 1 - [ ^ 



(4.24) 
(4.25) 

(4.26) 



For the functions (|4.19|) , l|4.20() we get here 
f(u) = u + 7700 

/'(«) = 1 - 



m 



(fJ-3 ~ u) 2 



(4.27) 
(4.28) 
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Inserting the solution k from 14.25|l . 14.26|l in Ji(k) we 
get 

/(0)= Mk)\ stat , (4.29) 
/'(0) = 1 - fe| > k\ > 0. (4.30) 

Clearly these arguments work similarly, if only one of 
the n a is equal to zero or all three u a are zero. In all 
cases H4.29fl holds for the exceptional points with |fc| < 
1, which can exist only if f(u) has no pole at u = 0. 
Since H4.29[) involves only " scalar" quantities, it holds in 
any basis. 

The case of exceptional solutions for |fc| — 1 can be 
treated in an analogous way. An exceptional solution 
of H4.10fl with u = fi a (a = 1, 2, 3) can only exist if the 
corresponding rj a = 0. Then f(u) has no pole for u = fi a 
and the exceptional solutions of (|4.1U|) fulfil 

fe = fe||+fe_L, (4.31) 

with 

fc|| =-(E-u)- 1 V \ u= ^, (E-^ a )k ± =0 

and 

/( Mo )= J 4 (fe)| stat , (4.32) 
f'(u a ) = 1 - fcf = fei > 0. (4.33) 

Note that if a solution is possible, fcj_ may be any linear 
combination of the eigenvectors to the eigenvalue u a of 
E, where the norm is given by |fej_| = y/ f'(/i a ). Thus we 
see that the function f(u) is very useful for discussing the 
stability of the THDM potential. What we have shown 
so far can be formulated as follows. 

Consider the functions f(u) and f'(u). Denote by /, 

/ = {ui,... s «n}, (4.34) 

the following set of values of u. Include in / all u where 
f'(u) = 0. Add u = to I if /'(0) > 0. Consider then 
the eigenvalues u a (a = 1, 2, 3) of E. Add those fi a to / 
where f(u a ) is finite and f'(ii a ) > 0. We have n < 10. 
The values of the function J^k) at its stationary points 
are given by 

MV\ stat =f(ui) (4.35) 

with Ui G /. The potential is stable if f(ui) > for all 
Ui G /. Then the stability is given solely by the quartic 
terms in the potential. The potential is unstable if we 
have f{ui) < for at least one Ui G /. If we have f(ui) > 
for all Ui G / and f(ui) = for at least one Uj G I we 
have to consider in addition J2{k) in order to decide on 
the stability of the potential. 

We turn now to this latter case. We shall show that 
we have then to consider in addition the function 

9{u) ^Ho-^iE-u)- 1 ^ (4.36) 



for the stationary points of J&(k) with 

Mk)\ stat = f(u*) = 0. (4.37) 

We have for the vectors k satisfying l|4.37|l 

J 2 {k)^g{u i ) (4.38) 

if Ui ^ fi a , that is Ui is not an eigenvalue of E. If Ui 
is an eigenvalue of E, that is Uj = u a G /, and ei{ui) 
(1 = 1,..., N) are the N < 3 eigenvectors to Ui, then we 
have 

inf J 2 (fe) = g(ui) - \£ ± (!h)\ Vf'M, (4.39) 

k 

where the infimum is taken over all exceptional solu- 
tions k to Ui and 

^W : =En^ e 'W- ( 4 - 4 °) 

We summarise our findings in a theorem. 

Theorem 1. The most general potential of the two- 
Higgs- doublet model has the form l|3.12|l . Its stability is 
decided in the following way. If the potential has only the 
quadratic termVi, it is stable for £o > marginally sta- 
ble for £o = |£| and unstable for £o < Suppose now 
that V4 ^ 0. We construct then the functions f(u) (|4.15() . 
f'(u) $t]]$ and g(u) (|Q6l . and the set I of (at 

most 10) u values. 

(i) If fiui) > for all Ui £ I the potential is stable in 
the strong sense l|4.6l) . 

(ii) If f(ui) < for at least one Ui G / the potential is 
unstable. 

(Hi) If f(ui) > for all Ui G / and f(ui) = for 
at least one Ui G / we consider also the function 
g(u) H4.36fl . The potential is stable in the weak 
sense (|4.5() if for all Ui G I where f(ui) = the 
following holds (see (|4"3St to (|PU|) ): 

g(ui)>0if Ul ^ [i a , (4.41) 
g(u.i) - \£ x ( Ui )\ y/f'( Ui ) > ifm = a a . (4.42) 

// in some or all of these cases we have = in- 
stead of > we have marginal stability l|4.4|l . // in 
at least one case we have < instead of > the 
potential is unstable. 

Our theorem gives a complete characterisation of the 
stability properties of the general THDM potential. In 
the following subsection we apply the theorem to assert 
the strong stability condition l|4.6|l for a specific potential. 
An application for the weaker stability condition l|4.5|) is 
given in Subsection 18. II 
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4.1. Stability for THDM of Gunion et al. 

We consider the THDM of with the Higgs potential 

V(fi,(p 2 ) = \i(<p\<pi - vf) 2 + X 2 (flf2 - vl) 2 
+ A 3 



+ A 4 
+ A 5 
+ A 6 
+ A 7 



which contains nine real parameters if we do not count 
the constant. This potential breaks the discrete symme- 
try 



•■P2 



<P2 



(4.44) 



only softly, i.e. by V 2 terms, thus suppressing large flavour 
changing neutral currents. For various restrictions on the 
THDM by symmetries see for instance 35] . Dropping the 
constant term, we put the potential into the form (|3.12|) 
using the relations (|3.11[1 . Then, 



??oo = t(Ai + A 2 +4A 3 + A 4 ), 






Ai — A2 



/2(A 5 -A 4 ) A 7 
E = - A 7 2(A 6 - A 4 ) 
8 \ 



From (|4.15|) and (|4.16f) we obtain 






2(Ai +A 2 - A 4 ) 



(4.45) 



f(u) = u+-{\ 1 + \ 2 + 4A 3 + A 4 ) 
(Ai - A 2 ) 2 



4(Ai + A 2 - A 4 - 4m) ' 



(4.46) 



(A1-A2)' 



/'(«) = 1 - 

(Ai + A 2 - A 4 - AuY 

We introduce the abbreviation 
1 



« : = g ( Afi + Ae ~ V ( As ~ Ae ) 2 + A ? 



Applying theorem^to the functions /(it) and f'(u), we 
find the strong stability assertion by V4, see (|4.6[) . to be 
equivalent to the simultaneous conditions 



Ai + A 3 > 0, 



A 2 + A 3 > 0, 



A 4 , k > -2A 3 - 2V(Ai+A 3 )(A 2 + A 3 ). 



(4.48) 



Re((/?|y 2 ) -v^cosi 

hn.(tp\(p 2 ) - ui« 2 sin^ 

Re((p\(p 2 ) -Ui^cos^ (lm((p\ip 2 ) -viv 2 sin£j , 

(4.43) 



In particular, if A 4 , A 2 , A 3 , A 4 , k > these inequalities are 
fulfilled. They can then be rewritten as: 

Ai,A 2 ,A 3 ,A 4 > 0, 4A 5 A 6 >A 7 ! . (4.49) 

That means, if (|4.49|l holds, the conditions l|4.48|l are ful- 
filled and the potential is stable. For the case A7 = 
we can replace k by min(A5,Ae) in the stability condi- 
tions H4.48fl . which arc then in particular fulfilled if A,j > 
for i = 1, . . . , 6. 

The potential in [9j is even more specific since it is in- 
variant under 14.44fl . Inserting their potential parameters 
in J4.48JI we reproduce the result of 9], their equation (2). 



5. LOCATION OF STATIONARY POINTS 

After our stability analysis in the preceding section 
we now determine the location of the stationary points 
of the potential, since among these points there are the 
local and global minima. To this end we define 

*-(£)■ *-(!)• MT*)- (5 ' 1) 

In this notation the potential l|3.12|l reads 

V = K T £ + K T EK 
and is defined on the domain 

K T gK > 0, K > 0, 



with 
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1 

-1 



(5.2) 
(5.3) 
(5.4) 



For the discussion of the stationary points of V, we 
distinguish the allowed cases K = 0, Kq > \K\ and 
K = \K\ > 0. 

The trivial configuration K = is a stationary point 
of the potential with V = 0, as a direct consequence of 
the definitions. 

The stationary points of V in the inner part of the 
domain, Kq > \K\, are given by 

1 - -t - 

EK = --£ with K gK > and K > 0. (5.5) 

For det E ^ we obtain the unique solution 



(5.6) 



(4.47) provided that 



-T - 



£ E gE £ > and K > 0, (5.7) 

and no solution if 15. 7|) does not hold. The Hessian ma- 
trix 

d 2 



dKi dKj 



V)=2E, where i,j = 0...3, (5.8) 
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determines whether (|5.6|l is a local minimum, a local 
maximum or a saddle. In the case det E = we may 
have exceptional solutions of Ij5.5(l . In the regular case 
as well as in the exceptional cases, the existence of a 
solution of (|5.5(l along with the corresponding values of 
the potential are not affected by the transformation of 
parameters (I3.18|) . 

The stationary points of V on the domain boundary 
Kq = \K | > are stationary points of the function 



F(K,u) :=V-uK gK, 



(5.9) 



where u is a Lagrange multiplier. The relevant stationary 
points of F are given by 

1 _ -t - 

{E - ug)K = --£ with K gK = and K Q > 0. 

(5.10) 

For regular values of u with det(-E — ug) / we obtain 
k{u) = - 1 ~(E-u~g)- 1 i (5-11) 

The Lagrange multiplier is determined from the con- 
straints in (|5.10|l by inserting (|5.11(l : 

g T (E-ug)~ 1 g(E-ug)~ 1 £ = and K > 0. (5.12) 

There may be up to 4 values u = fi a with a = 1, ... ,4 
for which det(E — ug) = 0. Depending on the potential 
some or all of them may lead to exceptional solutions 
of (|5.10|) . Note that for the regular as well as for the ex- 
ceptional cases, the Lagrange multipliers u and the value 
of the potential belonging to solutions (u, K) of (|5.10|l 
are, similar to above, invariant under the transforma- 
tions ipTT^i . 

For any stationary point the potential is given by 

V\ stat = ^K T £ = -K T EK. (5.13) 

Suppose now that the weak stability condition (|4.5() 
holds. Then (|5.13(1 gives for non-trivial stationary points 
where K ^ 0: 

V\ stat < 0, (5.14) 

since the cases V4 < and V4 — V2 = are excluded by 
the stability condition. 

Similarly to the stability analysis in Sect. 01 we can use 
a unified description for the regular stationary points of V 
with K > for both \K\ < K and \K\ = K defining 
the function 

f(u):=F(K(u),u), (5.15) 
where K(u) is the solution i|5.11|l . It follows: 

(5.16) 
(5.17) 



f(u) = -\i T {E-ug) l £ 
f{n) = -\( T {E~u~g)- 1 ~g{E- 



ug) £. 



Denoting the first component of K{u) as Kq(u) we sum- 
marise as follows. 



Theorem 2. The stationary points of the potential are 
given by 

(la) K = K{0) iff'(0) < 0, K Q {0) > and detE ^ 0, 

(lb) solutions K of (|5.5|l if det E = 0. 

(II a) K = K{u) for u with det(£ - ug) ^ 0, /'(it) = 
and Kq(u) > 0, 

(lib) solutions K of (|5.10|) for u with det(E — ug) = 0, 

(III) K = 0. 

In many cases, for instance if all values /2i, . . . ,/X4 are 
different, we can diagonalise the in general non-hermitian 
matrix gE in the following way: 



gE = Y,~^a- 



(5.18) 



Here the P a are quasi-projectors constructed from the 
normalised right and left eigenvectors Xa , Xa 

oigE. We 

have then gExa = fiaXa, XagE = Xap-a, XaXb = Sab 
and can impose as additional normalisation condition 
xlxa = 1- The P a are given by 



(5.19) 



and satisfy 



trP Q = l, P a P h = Fa fOTa = 6 ' (5.20) 

|0 for a ^ b, ( ' 

where a, b G {1, . . . , 4}. In terms of the P a (|5.16|) 
and H5.17fl read 



f n 1 

fm = — A 2^— ) 

4 ^ Ma - U 
a— 1 

4 ~T ~ 



(5.21) 
(5.22) 



Of course, f(u) (|5.16() is always a meromorphic function 
of u but in general poles of higher order than one may 
also occur. 



6. CRITERIA FOR ELECTROWEAK 
SYMMETRY BREAKING 

The global minimum will be among the stationary 
points discussed in the previous section. Here we dis- 
cuss the spontaneous symmetry breaking features of the 
possible classes of minima and give criteria to ensure a 
global minimum with the required electroweak symmetry 
breaking SU(2) L x U(l) Y _ -> U(l) em . 

A global minimum at K = means vanishing fields for 
the vacuum. In this case, no symmetry is spontaneously 
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broken. If the global minimum lies at K ^ 0, the full 
gauge group or a subgroup is broken. We denote the 
vacuum expectation values, i.e. the fields at the global 
minimum of the potential V, by 



(6.1) 



with i = 1,2. In general the vf 7 v® are complex numbers. 
To exhibit the consequences of electromagnetic gauge in- 
variance we consider the matrix l|3.8|) at the global min- 
imum, K | min . 

If the global minimum of V occurs with Kq > \K\, it 
follows detif| min > 0, see Sect. Since we have 



detK\ 



,,o „,+ l 



the vectors 



(6.2) 



(6.3) 



are linearly independent. Then there is no transforma- 
tion l|3.13[l such that both and i>^~ become zero. This 
means that the full gauge group 5J7(2)l x U(l)y is bro- 
ken. We can also show this using the methods explained 
in Appendix El see ljA23|) . 

In the case that the global minimum of V features 
Kq = \K\ > 0, the rank of the matrix K_\ m i n is 1 and the 
vectors <|6.3[1 are linearly dependent. After performing a 
SU(2)l x U(l)y transformation we achieve 



= 0, 



^0, 



(6.4) 
(6.5) 



and identify the unbroken U(l) gauge group with the 
electromagnetic one. By a transformation i|3.13|l . namely 



cos (3 sin j3 e ^ 
— sin p cos 



(6.6) 



with P fulfilling | ^j 1 1 sin /3 = [t^cos/?, we can arrange 
that 



»5° 



o. 



v /V2 




v > 



(6.7) 
(6.8) 



for the vacuum expectation values in the new basis. This 
can also be derived from the results in Appendix ^ 
see (|A16j) et seq. In (I6.8|) vo is the usual Higgs-field 
vacuum expectation value, vo ~ 246 GeV (see for in- 
stance [33|l. 

Now, we want to derive conditions for the parameters 
in the general potential l|3.12fl . which lead to the required 
EWSB by a global minimum with Kq = \K\ > 0. In the 
following, we assume the potential to be stable. If we 



consider parameters fulfilling £0 > |£| this immediately 
implies J2(fe) > and hence from (|4.5|) V > for all 
K 7^ 0. Therefore for these parameters the global min- 
imum is at K = 0. Thus we arrive at the requirement 



Co < III- 



Here we obtain 



dV 



dK 



k fixed, 
K =0 



(6.9) 



(6.10) 



for some k, i.e. the global minimum of V lies at K 7^ 0. 

Addressing the non-trivial cases, suppose that the two 
points 



P 



(6.11) 



withpo ^ |p| and go ^ \q\ ar< 3 stationary points of V, that 
is each of them is either a solution of l|5.5|l . or, together 
with an appropriate Lagrange multiplier u p or u q for p 
or q, respectively, a solution of l|5.10|l . 

Firstly, consider po — \p\. From (|5.2I) and (|5.1UI) we 
have 



dV 



dK 



K fixed, 



Co + 2(Ep) = 2u p p . 



(6.12) 



If Up < 0, there are points K with Kq > po, K — p and 
lower potential in the neighbourhood of p, which there- 
fore can not be a minimum. We conclude that in a the- 
ory with the required EWSB the global minimum (which 
needs to be on the light cone) must have Lagrange multi- 
plier uq > 0. As we shall show in Sect.[7| the case uq = 
leads to zero mass for the physical charged Higgs field at 
tree level. This is unacceptable from a phenomenological 
point of view if we disregard the possibility of very large 
radiative corrections. Therefore, we find as condition for 
an acceptable theory 



u > 



(6.13) 



for the Lagrange multiplier uq of the global minimum. 
Secondly, for p — \p\ and q = \q\ we have from (|5.13|) 
and l(5~TU|> : 



1 



1 



V(p)-V(q) = -p^--rt 

= p T (u q g~- E)q- q T (u p g- E)p 

= (u q ~u p )p T gq. (6.14) 



Since p and q are vectors on the forward light cone, p T gq 
is always non-negative and zero only for p parallel to q. 
Furthermore, the case that two different p, q are parallel 
can not occur, since then l|6.14|l requires V(p) — V(q), 
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while (|5.13|l and (|5.14|) imply V(p) ^ V(q) for that case. 
Therefore we conclude 



V(p) < V(q). 



(6.15) 



Assuming po — \p\ and qo > \q\ we get from Ij5.13l) and 



V(p) - V(q) = -u p p T gq 



and 



V(p)-y(g) = (p- 



q) T E(p 



q)- 



(6.16) 



3.17) 



The first equation implies in particular, that a stationary 
point on the domain boundary with positive Lagrange 
multiplier will have a lower potential than any stationary 
point with Kq > K. From the second equation follows in 
this case that E has a negative eigenvalue. Since for the 
stationary point q in the interior of the light cone the 
Hessian matrix is 2E (see (15.811 s ). we see that q cannot 
be a local minimum. This result and the hierarchies of 
the stationary points derived above agree with [l4|. We 
summarise as follows. 



minimum with the 
symmetry breaking 
and absence of zero 



Theorem 3. A global 
spontaneous electroweak 
SU(2) L x U(1) Y -» U(l) em 
mass physical charged Higgses 

(I) requires £ < 

(II) is given and guaranteed by the stationary point of 
the classes (Ila) or (lib) of theorem with the 
largest Lagrange multiplier uq > 0. 

We remark that for two different stationary points in 
the inner part of the domain or with u = on its bound- 
ary, any linear combination of them with Kq > \K\ is a 
stationary point as well. These points therefore belong to 
one connected set of degenerate stationary points. Stabil- 
ity requires that this set contains points with Kq > \K\ 
and is bounded by points with Kq = \K\. If interpreted 
geometrically, this degenerate set is a line segment, el- 
lipsoidal area or volume. Together with the arguments 
above we find the following mutually exclusive possibil- 
ities for local minima, expressed in terms of the gauge 
invariant functions: one or multiple solutions with the 
required EWSB (Kq = \K\) or the aforementioned de- 
generate set of solutions (Kq > \K\) or one charge break- 
ing solution (Kq > \K\) or the trivial solution (K = 0) 



7. POTENTIAL AFTER ELECTROWEAK 
SYMMETRY BREAKING 

We assume a stable potential which leads to the desired 
symmetry breaking pattern as discussed in the previous 
sections and derive consequences for the resulting phys- 
ical fields in the following. We choose a unitary gauge 



and the basis for the scalar fields such that for the vac- 
uum expectation values relations (|6.7|l and (|6.8|l hold, 
and furthermore the fields satisfy 



Imp? (a;) = 0, 
Rep? (a;) > 0. 

We introduce as usual a shifted Higgs field 
p'(x) := V2 Reipl(x) - v . 

Then the two Higgs doublets are 



if!(x) 



1 



( 







<P2(x) = r i> { 



(7.1) 
(7.2) 
(7.3) 

(7.4) 
(7.5) 



V2 \y + p'(x) / 

In addition to p' there are two more neutral Higgs fields 

h! := V2 Re y%, h" := V2 Im <p% (7.6) 

and the charged fields 

H+:=<pt, H~~(H+y. (7.7) 

It is convenient to decompose K according to the power 
of the physical fields they contain 



K = K {0} +K {1} +K {2} , 



with 




K {1} = v 



K 



(7.8) 



(7.9) 



2H_H+ + h' 2 +h"^ 



K 



{2} 



2p'ti 
2p'h" 



(7.10) 



vp' 2 -2H-H+ -h' 2 -h" 2 ' 



By Uq we denote again the Lagrange multiplier corre- 
sponding to the global minimum of V . From (15 . L Up we 
have 

M^^tio^o}-^". (7.11) 

From the explicit expressions (|7.9I) and l|7.10|l we further 
have 

4)^{orC 4}^{i} = °- ( 7 - 12 ) 



Using H7.8|) to l|7.12|l we obtain for the potential ijBT: 

V = V {0} + V {2] +V {3} +V {4} , (7.13) 

where Vs^\ are the terms of A; th order in the physical 
Higgs fields 

^{0} = (Co + 6) «g/4, (7.14) 

V {2} =K T {1} EK {1] +2uQK T { Q } ~gK {2} , (7.15) 



V {3} =2K {1} EK {2} , 

v {4} = kJ 0X ek 



{2}^^{2}- 



(7.16) 
(7.17) 
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The second order terms <|7.15[1 determine the masses of 
the physical Higgs fields: 



V {2} = -( P ',h',h")M 



2 

neutral 




m 2 H ±H + H~ 



(7.18) 



with 



M 



neutral 



-fo- 
-6 



6 



Vq V12 



-6 

Vn) v 2 r] 12 



2 M II 



V2 2 )j 

(7.19) 
(7.20) 



Note that the condition uq > corresponds to the pos- 
itivity of the charged Higgs mass squared at tree level. 
This result was already mentioned in Section |SJ Gener- 
ically the mass terms <|7.18[1 contain 7 real parameters. 
From l|7.19(l and (|7.20(1 we see that all 7 parameters are 
in general independent in this model. 



8. EXAMPLES 

Here we apply the general considerations of Sections 0] 
to0to specific models. 



8.1. MSSM Higgs potential 

In this subsection, we consider the MSSM Higgs poten- 
tial and reproduce the well-known results for its stability 
symmetry breaking and mass spectrum (see e.g. [24| and 
references therein), employing the method described in 
the previous sections. In the notation of ^jjj the MSSM 
Higgs potential is 



V = V D + V F + V S o{t 

with 

V D = \{g 2 +g l2 )(HlH 1 -HlH, 



.1) 



2 9 



H\H 2 



V f = \h\ 2 {H\h x + h\h 2 ), 
Koft = m 2 Hl + m 2 H2 h\h 2 - (ml H?eH 2 + h.c.) , 



where H\ and H 2 are Higgs doublets with weak hyper- 
charges y = —1/2 and y = +1/2, respectively, m 2 Hi , 

m 2 H2 E R, rn 2 6 (D and \[i\ 2 € Rq are parameters of di- 
mension mass squared. Substituting Hi and H 2 by dou- 
blets <f i, ip 2 with the same weak hypercharge y = +1/2 
according to l|3.5|) and using the relations l|3.11|l . we can 
put the potential in the form l|3.12|l . The parameters are 



for V4 = Vr> and 









-i( 








Co = |m| 2 + \( m Hi + m H 2 )> t = 



-Re(m§), 
Im(m§), 



'M 2 



(8.4) 

for V 2 = V F + V soft . 

We determine the stability of the potential by employ- 
ing theorem^] The functions f(u) I4.15|l and f'(u) (|4.1(j|) 
for the MSSM are 

f(u) = u+ ig 2 , (8.5) 
/'(«) = 1- (8.6) 
is given here by u = and the eigenvalues 



The set / 
of E (IQl . 



u\ = 0, u 2 



g 2 , u 3 



,'2 



(8.7) 



We find for the stationary points of J4 with ui — u\, U3 
the values J4(k)\ stat — f{ui) > but for those with 
u 2 the value J4(k)\ stat — f{u 2 ) = 0. Explicitely, the 
stationary points of J4 with u 2 are 

k = (cos cj>, sin 0, 0) T , (f> G R, with J 4 (fe) = 0. (8.8) 

They are known as the "D-flat" directions, since they 
have V D = 0. For the MSSM, they prevent the stability 
assertion by the quartic terms alone. For the stability to 
be guaranteed by V 2 > in these directions, theorem 
gives as condition, see <|4.42ll and (|4.36|) . the inequality 



sM - |€±M| v r FW) = Co -^Q + e 2 > 0. (8.9) 

Inserting (|8.4|) we get 

(8.10) 



\ m l\ < \vf + ^K, + m H 2 )- 



as the necessary and sufficient condition for the stability 
of the MSSM potential in the sense of (|4.5jl . 

For the global minimum to be non-trivial, criterion (I) 
of theorem |31 gives £0 < v£i + £2 + or equivalently 



1 



lA' + ~(™ Hl + m Z H2 ) < W \m 2 \ 2 + - (m 2 Hi - 



1 



'Hi 



(8.11) 

as a necessary and sufficient condition. We consider 
the acceptable global minimum candidates given by 
the classes (Ha) and (lib) of theorem [21 The condi- 
tions (|8.10|l and l|8.11|l prevent exceptional solutions. The 
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regular solutions are determined by the functions 



/(«) = 

/'(«) = 
2T (u) - 



iff 2 
8^ 



_if 

h l2 -u 



£o £l £2 



el 



(i5 2 
Co 

V - 

8^ 



,) 2 (- 



iff' 2 



w) 2 



(8.12) 
(8.13) 
(8.14) 



where we omitted the insertions H8.4(l for a compact no- 
tation. Employing again the conditions l|8.10|l and (|8.11l) 
we find the following. The function /'(it) always has two 
zeros and those zeros imply values of Kq(u) with oppo- 
site signs. The physical solution with Kq(u) > has the 
Lagrange multiplier 



!/() 



i l6|g 2 + y^-^-e 



2 V 2 



161 



(8.15) 



which is positive. Figure [2] shows the functions 
f'(u),Ko(u) for an example set of parameters (corre- 
sponding to the SPSla scenario |46j at the tree level). We 
conclude that l|8.10|l and l|8.11|) guarantee the existence of 
the stationary point K(uq), which fulfils criterion (II) of 
theorem [3] and therefore is the global minimum with the 
required EWSB pattern. Moreover, there are no other 
local minima. 

Note from (|8.3|) and (|8.4|) that it is always possible 



to choose a basis with £i = 



'3 h 



£2 = without af- 



fecting any other parameters. We further choose a gauge 
where (|6.4|) and (|6.5() hold, and perform the rotation (|6.6I) 
with 



tan/3 = 



/ gojM±vjO[Zllg 
Zq\z»\-V%-&-8& 



(8.16) 



into a basis of the form l|6.7|l and ^6.8(1 . which has the 
new parameters 



E' 



/-P2/S 


ml 




\-S2/3 


ml 


+ 




+ s 2/3 5 











g 2 




.17) 



(8.18) 



with the abbreviations g 2 :— g 2 + g 12 and S2/3 := sin 2/3 
etc. We insert the expressions into the formulae of Sec- 
tion and use 



\gv 



m 2 z 



\~gvo 
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FIG. 2: The global minimum determining functions f'(u) 
and K Q (u) for the MSSM, see lETTS! and lEHIt . with |^| 2 + 
m 2 Hd = 157486 GeV 2 , |^| 2 + m 2 Hu = -2541 GeV 2 , |m§| = 
15341 GeV 2 . The small boxes show the functions with en- 
hanced ordinate resolution in the region around the physically 
relevant zero of /'(«)■ 



within the Higgs sector of the MSSM. We obtain the 
mass squares 



to = 2w 2 (t7 22 + u Q ), m 2 H± = m 2 A0 + m\y 



.20) 



for the pseudo-scalar boson A := h" and the charged 
bosons H ± , which are already mass eigenstates. The 
non-diagonal part of the neutral mass matrix is 



Ml neutral, 
CP even 



4 rn% -hyfz ) (8 . 21) 

-2 s V m z m A° + s 2p m zJ 



in the basis (p,h'). Its diagonalisation leads to mass 
squares 



2 1 

m h°,H° = 2 ' '"' 



2 

A 



=F 



J(m\ +m%) 2 -4<%pm\ m%j (8.22 



(8.19) for the mass eig enstates h°,H°. They are obtained from 



with v = y / 2K (u ). According to (|7.19|) . the fact 
that £ 2 = V12 = implies tree-level CP conservation 



H° 
h° 



cos a' sin a' \ / p 



sin a cos a 



.23) 



15 



with the mixing angle a' determined by 



cos 2a' = 
sin 2a' = 



c 4 /3 m| 



S 4 /3 m| 



(8.24) 



Performing the rotation l|6.6[l on the complete doublets, 
as described above, leads to states p, h with simple cou- 
plings to the gauge bosons, e.g. vanishing ZZh' and 
WWh' couplings at tree level. However, usually the real 
parts of the neutral doublet components are excluded 
from that rotation. Applying the inverse of the rota- 
tion l|6.6l) to only the neutral components (p, h') gives 
\[2 (Re H\ , Re H 2 ) . The mass matrix in this basis is di- 
agonalised analogously to l|8.23[) . where a' is replaced by 
the mixing angle a with 



cos 2a = — cos 2(3 



sin 2a 



sin 2(3 



m 2 A o 


-m% 


m 2 H0 


-m 2 h0 






m A° 


+ m 2 z 


m 2 H o 


-m 2 



(8.25) 



which is the well-known result. 



8.2. Stationary points for THDM of Gunion et al. 

We continue the discussion of the potential from [l2l| . 
for which we derived the stability conditions in Section^ 
see (|4.43[) et seq. Note, that we consider the shifted po- 
tential according to V(K = 0) = 0. The parameters of 
V4 are given by 7700 , f] and E as in (|4.45|) , while we have 



£0 = -\\v\ - X 2 v 2 - 2X 3 (vl + v\ 
(-V1V2 (A 5 cos£ + 4f sin£)^ 
-V1V2 (A 6 sin£ + 4^ cos ^) 
-\iv\ + X 2 v 2 



\ 



(8.26) 



for the parameters of V 2 ■ Here, v\ , v 2 and £ denote 
the parameters of the potential (|4.43l) irrespective of 
their meaning for the vacuum expectation values. The 
function f'(u) constructed for this potential according 
to H5.17fl exhibits the zero 



M = ^A 4 , 



.27) 



whose associated solution K(p) is always a stationary 
point. It can be represented by the field configuration 



Pi la 





v 2 e*' 



(8.28) 



If all Xi > it is immediately obvious from the definition 
of the potential that this is the global minimum, which 
is furthermore non-trivial unless both vi,v 2 are zero. 



However, the quartic stability conditions l|4.48|l do not 
require the A^ to be positive. In the general case there 
can be more than one local minimum with the required 
EWSB, and the stationary point (|8.28l) is not necessarily 
the global minimum (we find parameters where l|8.28[) is 
only a saddle point, while another solution provides an 
admissible global minimum for a stable potential). The 
function f'(u) may have up to 5 additional zeros which 
can lead to further regular stationary points. Also excep- 
tional stationary points may occur for special parameter 
combinations. We do not find analytical expressions for 
the remaining zeros of f'(u) for the general case. Instead 
we apply the methods described in the previous sections 
in a semi-analytical way and determine the zeros of f'(u) 
numerically. We assure the stability conditions l|4.48|l to 
hold for the chosen parameter values. The stationary 
points with K = \K\ > and the largest Lagrange mul- 
tiplier u > are the required global minima. In order to 
classify also the other solutions, we compute m 2 H± and 
the eigenvalues of A / l^ cutral for each stationary point, as 
described in the previous sections for the global mini- 
mum. A solution for which all of these values are pos- 
itive is a local minimum. Mixed positive and negative 
values mean that the solution is a saddle point. Note 
that for a global minimum different from (|8.28|) . the po- 
tential parameters Vi, v 2 and £ lose their simple meaning. 
For example, we find that the global minimum acquires 
a non- vanishing, CP violating, phase for certain param- 
eters with £ = 0, A5 = A 6 7^ and A7 7^ 0. 

Figure [21 shows the potential at the stationary points 
with Kq = \K\ > for the particular parameter val- 
ues described in the caption. This example features tree- 
level CP conservation within the Higgs sector: £ = A7 = 
leads to £2 = and, by (|5.11l) and A 4 7^ A 6 , to K 2 = 
(i.e. trivial phases) at the global minimum. In the basis 
with (|6.7|) and (|6.8|l we then have £ 2 = V12 = 0; implying 
CP conservation by (|7.19|l . Note that even with the sim- 
ple parameter combinations chosen for Fig. |3| the struc- 
ture of the stationary points can be non-trivial. In the 
example Ai and A2 are equal and varied simultaneously. 
For the plotted range, where Ai is negative, the global 
minimum is a regular solution and differs from Q8.28[) . 
which is only a local minimum. For Ai = there are two 
exceptional and degenerate minima. The figure shows for 
positive Ai the expected behaviour, that l|8.28|l becomes 
the global minimum, but also that for < Ai < 0.0268 
a second regular local minimum exists. Two stationary 
points disappear for the plotted range above Ai > 0.0268 
because there the corresponding two zeros of f'(u) have 
a non- vanishing imaginary part. 



9. CONCLUSIONS 

We have analysed the scalar potential of the general 
THDM. In order to give an acceptable theory, this po- 
tential has to obey certain criteria. The potential should 
be stable, that is bounded from below, and lead to the 
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FIG. 3: The potential V of [l2T ]. shifted to V(K = 0) = 0, at all stationary points with Kq — \K\ > in dependence of Ai, 
where A 2 = Ai. The other parameters are A3 = 0.1, A 4 = 0.2, A 5 = A 6 = 0.4, Vi = 30 GeV, v 2 = 171 GeV, A 7 = 0, £ = 0. The 
lines represent regular stationary points, where the solid curve corresponds to the "obvious" solution with tan (3 — V2 jv\ and 
vo = \J2(y\ + «!), which is a local minimum for the chosen parameters. For Ai = A2 = there are two degenerate exceptional 
minima and regular solutions only for the saddle points. Depending on Ai = A2, the global minimum is given by that local 
minimum out of the two which has the lower potential. 



EWSB pattern observed in Nature. The conditions found 
for the stability of the potential and for EWSB are trans- 
parent and compact if the potential is written down in 
terms of gauge invariant field functions l|3.9|l . These 
conditions allow for a simple application to every spe- 
cific THDM. We illustrated our method in two examples, 
namely the MSSM potential as well as the THDM poten- 
tial introduced by Gunion et al. 12]. In the case of the 
MSSM we could easily reproduce all well-known results. 
For the potential of Gunion et al. we could clarify some 
interesting aspects of the model, see Figure |3 

We note that the method presented here may also be 
extendable to general Multi-Higgs-Doublet Models. A 
first step in this direction is done in Appendix [BJ Fur- 
ther, in a more detailed study it is mandatory to take 
quantum corrections to the Higgs potential into account. 
For the resulting effective Higgs potential the conditions 
for stability and for EWSB are then in general modified. 
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APPENDIX A: GAUGE ORBITS 

Here we give the proof that the gauge orbits of the 
Higgs fields of the THDM are parametrised by four 
vectors (Kq,K) satisfying i|3.19[l . or equivalently, by 
positive- semidcfinite matrices K_ (|3.8|l . (|3.9|l . Indeed, let 
us consider two Higgs doublets as in i|3.1fl . 



<p?(x), i = 1,2, a = +,0. 
We arrange these in a 2 x 2 matrix 

xi \ 1 at \\ ( <Pi( x ) <Pi(. x ) 



(Al) 



(A2) 



Then we have from (|3.8H 

K(x) =ct>(x)tf(x). (A3) 

The change of basis l|3.13|l means the transformation 

(f>(x) -> <f)'(x) = U<j){x). (A4) 

A gauge transformation from the SU(2)l x U(1)y gauge 
group means 

<0r) -> (x) = (U G (x)) a0 rf(x), (A5) 

where 

U G (x) g U(2). (A6) 
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Thus, under a gauge transformation the matrix <p(x) be- 
haves as 



(x)^<l>'(x)=cj>(x)U%(x). 



(A7) 



As we discussed in Section [3] any matrix K_{x) formed 
from the Higgs fields according to Ij3.8|l . which is equiv- 
alent to (|A3|) . must be positive semidefinite. Conversely, 
given any positive semidefinite matrix K_{x) we can di- 
agonalise it by a 2x2 unitary transformation W(x): 



K{ X )=W{x)^ K2 ° {x) )wHx), 
W\x)W{x) = 1. 
Since we have k.\{x) > and K2(x) > we can set 



(AS) 



««) = W(z) {ff> J^j (A9) 

and get 

K(x) =ct>(x)tf{x). (A10) 

With this we have proven the following. 

• For any positive semidefinite matrix Kjx) there are 
Higgs fields satisfying (|A3|) respectively (|3.8|) . 

Now suppose that we have a given positive semidefinite 
matrix K_{x) and two Higgs-field matrices <p{x) , 4>' {x) , 
both satisfying l|A3(l . 



K(x) = (f>(x)^(x) = (f>'(x)(f> n (x). (All) 

We want to show that 4>'(x) and 4>{x) are then related by 
a gauge transformation (|A7() . We consider three cases. 

(i) K(x) = 0. Then 4>{x) = cf>'(x) = and JSZJ is 
trivially fulfilled. 

(ii) K_(x) > 0, that is K_(x) is positive definite. Then 
detK{x) = \detct)(x)\ 2 = |det cf>' {x)\ 2 > (A12) 
and both <f)(x) and <^'(a;) have an inverse. We set 

(j>-\x)<f>'{x)=U%(x) (A13) 
and find from (|A11I) 

U G (x)U G (x) = 1, (A14) 

that is U G (x) £ U(2), and 

tf(x) = 4>{x)lf£{x). (A15) 

Thus (/''(a;) and 4>(x) satisfy (|A7|) . they are related 
by a gauge transformation. 



(iii) K_ (x) has rank 1, that is the eigenvalues are 

kiO) >0, k 2 (x) = 0. (A16) 

With the matrix W(x) diagonalising K_(x) 
(see (lA8fl ) we have then from IjAllll 

Kl n ] °r) = (wHx)cb(x)) (W\x)cb(x)y 



(wHx)(f>'(x)) (Wi{x)<t>'(x)) j 



From this we see that 

W\x )( f>(x) = (*t(*) x8(«)) , 

= ( x ' 1+ (x) x ' l0 (a;) ) , 

where 

XiOOxiO) = Xi{x)x'i(x) = Ki(x), 



(A17) 



(A18) 



(A19) 



with the composition of the vectors Xi( x )> Xi( x ) de- 
fined as for the vectors <pi (x) . Therefore we can find 
a matrix Uq(x) £ U(2) such that 



x'i a (x) = (U G (x)) a0X P(x), 

which implies 

W(x)cf>'(x) = W(x)(f>(x)U£(x), 
4>'{x) = 4>{x)Ul{x). 



(A20) 



(A21) 



That is 4>'(x) and (j>(x) are related by a gauge trans- 
formation. 

With this we have completed the proof of the following 
statement. 

Theorem 4. Any two Higgs-doublet fields giving the 
same matrix K(x) (|3.8|) . respectively (|A3p . are related by 
a gauge transformation. The space of gauge orbits can be 
parametrised by four-vectors (Kq,K) lying on and inside 
the forward light cone, see (|3.19() . 

We close this appendix with a supplementary note to 
Section|H|on the breaking of the SU(2)l x U(1)y gauge 
group. The matrix (|A2|) of vacuum expectation values is 



,+ ,,o 



(A22) 



If the global minimum of the potential occurs with 
Kq > \K | the corresponding matrix K_ has rank 2. Then 
from (|A3|) and (|A12|) we see that also 4> vac has rank 2. 
Invariance of tfi vac under a gauge transformation (|A7|) 



t>vac Uq 



(A23) 



is then only possible for Uq — 1. That is how we see 
with the methods of this appendix that in this case the 
full gauge group SU(2)l x U(1)y is broken. 
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APPENDIX B: THE CASE OF TV DOUBLETS 

In this appendix we generalise the methods of Sectional 
and Appendix lAlto the case of n > 2 Higgs doublets. We 
consider n complex Higgs-doublet fields 



tpi(x) 



1, . . . , n. 



(Bl) 



All doublets are supposed to have the same weak hyper- 
charge y = +1/2. In analogy to i|3.8[l we introduce the 
matrix 



(B2) 



which is now a n x n matrix. The aim is to discuss the 
properties of K(x) . For this we introduce the nxn ma- 
trix 4>(x) (compare HA2jl ) 



4>(x) 



(ip+{x) ifHx) 
<pt{x) ${x) 



0\ 


0/ 



(B3) 



\ V t(x) <p° n (x) 
It is easy to see that we have 

K(x) = 4>(x)^(x). (B4) 
A change of basis among the doublets means 

<j>{x)^4> l {x) = U^{x) (B5) 
with a constant matrix U £ U(n), 

U ] U=l n . (B6) 
A gauge transformation from SU(2)l x U(1)y means 

4>{x)^cj>'{x) = ( j>{x)Ul{x), (B7) 
where Ug(x) is block-diagonal 



U G (x) := 



fU G {x) 


MM 


v o 


l„-2/ 



(B8) 



U G {x) £ U{2) and thus U G {x) £ U{n). We have then 
from JET 



d°>(x) = (U G (x)) a0 cpf(x), 



1. 



(B9) 



From l|B3|) and (|B4|I we see that the matrix K(x) has the 
following properties 



• K.{ x ) i s positive semidefinite, 

• K.( x ) has rank < 2. 



That is, K_{x) has at most two eigenvalues 
Ki(x), k 2 (x) > and the remaining eigenvalues 
K3(x), . . . , K n (x) must be zero. The rank condition 



can be seen as follows. We denote by ip + (x), ip°(x) the 
first two column vectors of 4>{x). Then we have 



4>{x) = (^ + (x),^°(x),0,...,0j, (BIO) 
K{x) = L+*(x) 4,+ {x) + <pl*(x) i>°(x),..., 

^*{x)i>+{x) + ^*{x)^\x)\. (Bll) 



That is, at most two column vectors of K_(x) are linearly 
independent. 

Suppose now that we have a given positive semidefinite 
matrix K_(x) of rank < 2. Then we can diagonalise K(x) 
and represent it as 



K{x) = W{x) 



( Kx{x) 
K 2 (x) 



V 



W\x) (B12) 



with W(x) £ U(n) and fti(x) > 0, n 2 {x) > 0. We set now 



(x) = W(x) 



fy/^{x) 

\Zk 2 (x) 



V 







(B13) 



0; 



and see easily that (f>(x) is of the form l|B3|l and satis- 
fies <|B4(1 . Thus to any positive semidefinite matrix K(x) 
of rank < 2 there is at least one field configuration of 
the n Higgs doublets such that l|B2|l holds. 

Suppose now that we have two field configurations, 
that is two matrices <fi(x) and 4>'( x ) of the type 1|B3J| such 
that 

K{x) =(f>(x)tf{x) =ct>'(x)(t>' i (x). (B14) 
We can diagonalise K_(x) as in 1|B12(1 and get 



( Kl{x) 

n 2 (x) 



V 







(W\x) <t>{x)) (W^ {x) (t){x))^ 

(W\x) (f/(x)) (WHx) cj>'(x)) ] . 

(B15) 



From this we see that we must have 



W\x)4){x) 



( xiix) 


x2(*) 





xt( x ) 




V 





'x[ + ( x ) X[°(x) 
W^x)c(>'(x)= I X ' 2 + ( x ) x' 2 °( x ) 



(B16) 



(B17) 
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where 

xl(x)x2{x) =X2( X )X2( X ) = K z(x), (B18) 
Xi(x) X2(x) = Xi{x)x' 2 (x) = 0. 

From this we conclude that we can find a ma- 
trix U G (x) e 17(2) such that 

X'i a = {U G {x)) ap xf (x), i = 1,2. (B19) 

Inserting this Ug{x) into l|B8(l we get 

W\x)(j)'{x) = W\x)tf>(x) U%{x) (B20) 

and since W(x) £ U(n) 

4>'{x) = 4>(x)U%(x). (B21) 



That is 0'(x) and cj>(x) are related by a gauge transfor- 
mation. We summarise our findings in a theorem. 

Theorem 5. Forn Higgs- doublet fields of the same weak 
hypercharge y = +1/2 £/ie matrix K(x) = up^(x)(pi(x)\ 
is a positive semidefinite nxn matrix of rank < 2. For 
any positive semidefinite nxn matrix K(x) of rank < 2 
there are Higgs fields such that QBfy) holds. Any two field 
configurations giving the same matrix K(x) are related 
by a SU(2)l x U(1)y gauge transformation. The matri- 
ces K(x) form, therefore, the space of the gauge orbits of 
the n Higgs- doublet fields. 

As an example we consider three Higgs doublets. The 
space of gauge orbits is then given by all positive semidef- 
inite 3x3 matrices Kfx) with det K_{x) = 0. 
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